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ABSTRACT 

This  progress  report  outlines  our  research  efforts  on  modeling,  analyses  and  optimization  of 
multi-layered  printed  antennas  for  high  gain  applications.  Our  research  during  this  interim 
period  has  resulted  in  progress  in  the  following  four  areas.  1)  Performance  analysis  of  a  linear 
array  ofYagi-like  printed  sub-array  antennas.  The  effects  of  mutual  couplings  on  gain  and  pattern 
degradation  of  previously  optimized  Yagi-like  structures  are  investigated  for  various  array 
configurations  that  one  may  encounter  in  a  practical  design.  2)  Development  of  an 
electromagnetic  optimization  engine  based  on  Method  of  Moments  and  Evolutionary 
Programming  for  design  of  multi-layer  printed  antennas.  During  this  interim  period,  this  code 
has  been  extended  to  include  printed  patch  antennas  of  arbitrary  shapes.  In  addition,  the  use  of 
various  mutation  operators  to  speed-up  the  optimization  process  has  been  investigated  in  details 
where  it  has  been  found  that  a  Cauchy  mutation  operator  can  significantly  speed-up  the  optimal 
design  of  antenna  structures  particularly  when  the  corresponding  objective  functions  include  many 
local  optima.  3)  Design  of  circularly  polarized  (CP)  high  gain  multilayer  antenna  elements.  The 
design  feasibility  of  a  high  gain  CP  microstrip  patch  antenna  is  demonstrated  for  MS  AT 
application  where  a  gain  of  more  than  11.6  dBi  with  an  axial  ratio  of  less  than  1.1  dB  is  obtained. 
4)  Theoretical  investigation  of  Gain  Enhancement  methods  for  the  Yagi-like  antennas  printed  in  a 
multi-layer  anisotropic  medium.  During  this  interim  period,  the  conditions  for  obtaining  high  gain 
together  with  high  radiation  efficiency  for  a  microstrip  dipole  antenna  embedded  in  a  multi-layer 
anisotropic  medium  have  been  derived.  In  particular,  it’s  shown  that  by  a  proper  selection  of 
anisotropy  ratios  and  thickness  of  the  layers,  it  is  possible  to  obtain  two  high  gain  beams  at  0E  and 
eH  in  E  and  H  planes  respectively.  In  addition,  the  physics  of  these  high  gain  fields  are  investigated 
in  terms  of  the  leaky-wave  radiation. 
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I.  INTRODUCTION 


This  progress  report  summarizes  our  research  efforts  on  modeling,  analyses  and  optimization 
of  multi-layered  printed  antennas  for  various  applications  in  communications  and  radar  that  may 
require  high  gain,  wide  bandwidth,  multi-band  and/or  circular  polarization.  The  report  covers  the 
time-period,  October  1998  -  September  1999.  The  technical  contributions  made  during  this  period 
include: 

i)  The  radiation  performance  of  Yagi-like  stacked  printed  sub-array  antennas  when  used  as 
elements  in  an  array  environment  have  been  analyzed.  In  particular,  the  effect  of  mutual  coupling  on 
gain  and  far-field  pattern  of  previously  optimized  Yagi-like  structures  are  investigated  for  various 
array  configurations  that  one  may  encounter  in  practice.  It’s  found  that  gain  degradation  due  to  the 
surface-waves  and  coupling  effects  can  be  substantial 

ii)  The  electromagnetic  optimization  engine,  which  its  development  was  reported  in  our  last  annual 
report  [1],  has  now  been  extended  to  include  multi-layer  printed  patch  antennas  of  arbitrary  shapes. 
This  numerical  engine  is  based  on  a  hybrid  combination  of  Method  of  Moments  (MOM)  and 
Evolutionary  Programming  (EP).  In  addition,  the  convergence  rates  of  two  different  mutation 
operators,  namely  Gaussian  and  Cauchy,  in  the  EP  optimization  process  have  been  investigated  in 
details.  It  has  been  found  that  the  Cauchy  mutation  operator  can  significantly  speed-up  the  global 
optimal  design  of  antenna  structures  whose  objective  functions  posses  many  local  optima  and  when 
a  few  or  no  constraints  are  put  on  the  design  parameters.  We  have  applied  this  optimization  engine 
to  optimize  the  gain  of  a  stacked  Yagi-like  array  of  narrow  strip  dipoles  as  well  as  Yagi-like  array  of 
rectangular  microstrip  patch  elements. 

ill)  In  order  to  demonstrate  the  design  feasibility  of  multi-layer  high  gain  microstrip  antennas  with 
Circular  polarization  (CP),  a  three-layer  rectangular  patch  antenna  is  designed  for  MSAT 
application.  A  MOM  analysis  of  this  structure  shows  a  gain  of  more  than  11.6  dBi  with  an  axial 
ratio  of  less  than  1.1  dB  at  an  operating  frequency  of  1.6  GHz.  In  addition,  the  mutual  coupling 
effects  on  gain  and  axial  ratio  of  these  antennas  in  an  array  environment  has  been  analyzed. 

iv)  As  a  first  step  toward  the  extension  of  the  Yagi-like  concept  to  antennas  printed  in  multi-layer 
anisotropic  media,  we  have  analyzed  the  conditions  for  obtaining  high  gain  together  with  high 
radiation  efficiency  for  a  microstrip  dipole  antenna  embedded  in  a  multi-layer  anisotropic  medium. 
It’s  shown  that  by  a  proper  selection  of  anisotrpy  ratios  and  thickness  of  the  layers,  it  is  possible  to 
obtain  two  high  gain  beams  at  0E  and  0H  in  E  and  H  planes  respectively.  Examples  of  various  3  and 
5  layer  configurations  have  been  studied.  In  addition,  the  physics  of  these  high  gain  fields  are 
investigated  in  terms  of  the  leaky-wave  radiation. 
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II.  PERFORMANCE  ANALYSIS  OF  A  LINEAR  ARRAY  OF  YAGI-LIKE  PRINTED 

SUB-ARRAY  ANTENNAS 

In  our  last  annual  technical  report  [1]  we  presented  the  analysis  of  a  microstrip  Yagi-like 
antenna  which  consists  of  a  reflector  (ground-plane),  a  driver  and  a  finite  number  of  embedded 
director  strip  elements  (Figure  1).  This  structure  gives  more  degrees  of  freedom  (i.e.,  the  dielectric 
constants  of  the  layers)  when  optimizing  its  design.  It  can  also  be  made  conformal  to  various 
shapes  and  surfaces.  In  addition,  unlike  the  gain-enhancement  techniques  reported  for  an  embedded 
microstrip  antenna  [2,3],  without  directors,  the  Yagi-like  structures  can  achieve  high  gains  without 
the  need  for  very  high  permittivity  dielectric  layers.  In  fact  we  have  previously  shown  that  it  is 
possible  to  optimize  a  three-layer  structures  with  dielectric  constant  of  er  =2  for  a  gain  as  high  as  10 
dBi  or  more.  The  Yagi  gain  enhancement  technique,  however,  results  in  narrow  gain  and  impedance 
bandwidths.  The  question,  therefore,  arises  on  how  these  Yagi-like  structures  perform  when  used  as 
radiating  elements  in  an  array  environment. 

M— _ ► 

dn 


d3 

d2 
dl 

Figure  1:  A  Yagi  stacked  microstrip  array  in  a  multi-layered  medium 

Let  us  consider  a  3  layer  microstrip  Yagi-like  antenna  with  er]  =  er2  =  er3  =2.  Table  I  shows 
the  design  parameters  for  three  optimized  cases  reported  in  the  last  annual  report  [1].  The  width  of 
each  strip  element  is  w  =  0.02  A0  where  A0  is  the  free-space  wavelength.  Case  1,  which  results  in 
more  than  13.5  dBi  of  gain,  was  optimized  with  a  constraint  on  the  input  reactance,  Xin  =0  (i.e., 
resonant  driver  element)  but  with  no  constraint  on  the  input  resistance,  Rin.  Solution  in  Case  2,  was 
obtained  subject  to  an  additional  constraint  of  Rin  >10  ohms,  while  in  Case  3,  in  addition  to  the 
above  constraints,  a  constraint  was  also  put  on  the  maximum  total  thickness:  D  =  d,  +  d2  +  d3  <  0.5 
A,d.  Gains  for  these  three  cases  as  functions  of  normalized  frequency  are  depicted  in  Figure  2,  while 
the  corresponding  E-plane  far-field  patterns  are  shown  in  Figure  3. 

We  now  examine  the  performance  of  these  three  cases  when  they  are  used  to  form  a  2- 
element  as  well  as  a  5-element  linear  array.  We  have  investigated  both  collinear  as  well  as  side-by- 
side  element  arrangements. 
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Table  1:  Previously  Optimized  3-layer  ‘Yagi-like’  stacked  microstrip  array 
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(dBi) 
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W 
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Case  1 

13.57 

0.240 

0.120 

0.360 

0.3180 

0.3295 

0.3889 

0.109 

93% 

Case  2 

11.10 

0.299 

0.252 

0.2955 

0.3005 

0.3327 

14.43 

96% 

Case  3 

9.88 

0.125 

0.092 

0.230 

0.2838 

0.2881 

0.3686 

20.32 

81  % 

Degree 


Frequency  Figure  3:  E-plane  far-field  patterns 

Figure  2:  Gain  versus  frequency 


A.  Collinear  Yagi-like  Arrays 

Figures  4  and  5  show  a  two  and  a  five  element  collinear  Yagi-like  array  with  an  inter-element 
spacing  of  do  =  0.5  X0,  respectively.  The  element  axis  and  array  axis  is  both  assumed  to  be  along 
the  x-axis.  Far-field  patterns  in  the  <|>  =  0  plane  are  plotted  when  the  radiating  elements  correspond 
to  one  of  the  three  optimized  cases  in  Table  1.  In  order  to  show  the  effect  of  mutual  couplings,  the 
far-fields  computed  when  the  couplings  are  neglected  (i.e.,  element  pattern  multiply  by  array  factor) 
are  included  for  comparison.  As  can  be  seen  the  coupling  has  minimal  effects  on  the  pattern  for  the 
2-element  array  except  perhaps  for  Case  1  which  as  was  shown  in  Figure  2  is  very  narrow  band  and 
therefore  more  susceptible  to  mutual  coupling  and  surface-wave  effects.  The  effects  of  mutual 
coupling  are  more  pronounced  for  the  5-element  array  in  Figure  5.  The  pattern  is  greatly  influenced 
when  the  radiating  elements  are  those  of  Casel  and  the  least  affected  by  the  coupling  effects  when 
the  elements  are  those  of  Case  3  which  has  a  relatively  broad  band  frequency  response.  One  can 
particularly  notice  the  large  increases  in  the  side-lobe-levels.  The  pattern  degradation  for  Case  1  is 
partly  an  indication  that  the  surface-wave  reduction  mechanism,  resulted  from  optimization  of  the 
lengths  of  the  elements  and  the  thickness  of  the  dielectric  layers,  has  been  disturbed  by  the  mutual 
coupling  effects. 
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Figure  4:  Patterns  in  the  <|>  =  0  plane  for  a  2-element  collinear  array  of  Yagi-like  antennas 
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Figure  5:  Patterns  in  the  <|>  =  0  plane  for  a  5-element  collinear  array  of  Yagi-like  antennas 


B.  Side-by-side  Yagi-like  Arrays 

Figures  6,  7  and  8  show  the  results  for  two  and  five  element  Yagi-like  arrays  when  the 
radiating  elements  are  positioned  in  a  side-by-side  arrangement.  The  element  axis  is  x-axis  and  the 
array  axis  is  assumed  to  be  along  the  y-axis.  For  the  five  element  array,  patterns  in  both  the  <|>  =  0 
and  the  <J>  =  90°  planes  are  plotted.  Again,  the  effect  of  mutual  coupling  is  more  pronounced  for 
Case  1  particularly  in  the  §  =  0°  plane  where  the  pattern  without  the  coupling  effects  is  basically  that 
of  a  single  Yagi-like  element  shown  in  Figure  3.  The  patterns  in  the  (j)  =  0  plane  for  Case  2  and 
Case  3  are  relatively  unaffected  by  the  coupling  effects  because  of  their  broader  radiating  elements 
frequency  responses, 

Finally,  gain  and  radiation  efficiency  (i.e.,  ratio  of  the  radiated  power  to  the  total  of  the 
radiated  plus  the  surface-wave  powers)  are  tabulated  in  Tables  2  through  4  for  the  three  Yagi-like 
element  cases.  Arraying  does  not  significantly  increase  the  gain  for  Case  1  because  of  the  strong 
coupling  effects.  Without  coupling  effects  one  would  expect  an  increase  of  approximately  3  and  7 
dBi ,  over  the  gain  of  an  isolated  element,  for  a  2-element  and  a  5-element  array,  respectively.  The 
results  in  these  tables  together  with  the  above  far-field  patterns  clearly  show  that  in  order  to  fully 
take  advantage  of  the  Yagi-like  gain  enhancement  technique,  one  should  globally  optimize  gain  of 
the  whole  array,  in  the  presence  of  the  mutual  coupling  effects,  rather  than  an  element  in  isolation. 


Table  2:  Case  1 


Single 

element 

2  elements 
collinear 

2  elements  side 
by  side 

5  elements 
collinear 

5  elements 
side  by  side 

Gain  (dBi) 

13.57 

14.41 

13.70 

16.07 

15.29 

Efficiency 

92.7% 

96.5% 

89.8% 

94.6% 

88.6% 

Gain  improvement 
over  single  element  (dBi) 

1.33 

1.34 

3.83 

4.42 

Table  3:  Case  2 


Single 

element 

2  elements 
collinear 

2  elements  side 
by  side 

5  elements 
collinear 

5  elements 
side  by  side 

Gain  (dBi) 

11.10 

12.43 

12.44 

14.93 

15.52 

Efficiency 

96.7% 

97.3% 

94.8% 

97.3% 

94.1% 

Gain  improvement 
over  single  element  (dBi) 

- 

1.33 

1.34 

3.83 

4.42 

Table  4:  Case  3 


Single 

element 

2  elements 
collinear 

2  elements  side 
by  side 

5  elements 
collinear 

5  elements 
side  by  side 

Gain  (dBi) 

9.88 

11.31 

11.55 

14.32 

14.86 

Efficiency 

81% 

77.1% 

83.2% 

76.9% 

84.9% 

Gain  improvement 
over  single  element  (dBi) 

~ 

1.43 

1.67 

4.44 

4.98 
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Figure  6:  Patterns  in  the  (|>  =  90°  plane  for  a  2-element  side-by-side  array  of  Yagi-like  antennas 
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Figure  7:  Patterns  in  the  $  =  90°  plane  for  a  5-element  side-by-side  array  of  Yagi-like  antennas 
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HI.  AN  ELECTROMAGNETIC  OTIMIZATION  ENGINE  FOR  MULTI-LAYER 
PRINTED  ANTENNA  STRUCTURES 

Designs  of  low  profile,  planar  and  conformal  antennas  in  a  multi-layered  substrate 
environment  are  rapidly  becoming  of  increasing  importance  in  many  radar  and  communication 
applications.  Examples  include,  design  of  stacked  patch  antennas  for  broadband  or  multiple- 
frequency,  dual-polarized  and  high  gain  applications.  To  optimize  and  design  these  antenna 
structures,  an  efficient  numerical  engine  for  electromagnetic  modeling  together  with  multi-objective  / 
multi-parameter  global  optimization  capabilities  is  needed.  There  is  no  treatment  of  a  general 
multiply  stacked  multiple  patches  (MSMP)  in  open  literature  though  theoretical  and  experimental 
analyses  of  stacked  patches  in  a  two-layer  microstrip  configuration  or  multiple  patches  on  a  single¬ 
layer  microstrip  substrate  are  available.  On  the  other  hand,  several  proprietary  software  packages 
such  as  Ensemble  [4]  and  IE3D  [5]  can  be  utilized  for  analysis  but  not  for  optimal  design. 

In  this  work,  we  have  developed  an  optimization  and  design  engine  which  uses  a  hybrid 
combination  of  the  method  of  moments  (MOM)  and  evolutionary  programming  (EP)  to  globally 
optimize  a  given  MSMP  antenna  structure. 

A.  Implementation  of  the  Method  of  Moments 

We  have  used  the  mixed-potential  integral  equation  (MPIE)[6,7]  to  model  planar  microstrip 
structures  of  arbitrary  shape  in  a  multi-layered  dielectric  medium  (Figure  9).  For  P(m)  patches  on 
mth  layer  as  shown  in  Fig.  9,  there  are  P(m)  MPIE’s  in  the  following  form  (p=  1, 2, . . .,  P(m )): 


VG, 


mp,nq 


(r  |  r  )V  •  J^fr  ) 


mpJ»p(r)  =  -*>>(*•) 


(1) 


Z 


y  * 


Figure  9:  A  multi-layered  microstrip  antenna  structure 
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In  (1),  and  Zmp  are  the  electric  surface  current  and  the  surface  impedance  of  pth  patch  on  mth 
layer,  respectively.  is  the  excited  tangential  electric  field  on  that  patch.  For  the  present  problem, 
it  is  assumed  that  an  ideal  voltage  gap  generator  is  connected  to  the  feed  at  the  truncated  end. 
G%P,nq  and  GEmpnq  are  the  scalar  Green’s  functions  of  electric  and  magnetic  types,  evaluated  at  pth 

patch  on  mth  layer  due  to  a  source  at  qth  patch  on  nth  layer.  These  functions,  which  are  of  the 
Sommerfeld  type  and  exhibit  a  translational  invariance  along  x  and  y,  can  be  pre-computed  for  a 
range  of  lateral  distance.  Our  implementation  of  the  MPIE  utilizes  a  novel  technique  for  a  semi- 
analytical  evaluation  of  these  Green’s  functions  and  is,  therefore,  well  suited  for  a  medium  with 
many  dielectric  layers  [1,8].  In  total,  there  are  P(l)+  P(2)+-  •  •+  P(N)  equations  like  (1). 

To  solve  the  equations  for  the  unknown  surface  current,  moment  method  with  Galerkin’s 
testing  procedure  has  been  used.  The  process  begins  with  the  following  expansion  of  the  x-  and  y- 
current  into  rooftop  basis  functions  in  x-  and  y-direction, 

U(nq)  V  (nq  ) 

ws;  (*'-*„,)  2X  <V  +v)^4  (2) 

u- 1  V=1 

where  Inq(u)  and  Inq(U+v)  are  the  coefficients  of  the  basis  functions  K  and  Bnq  5  respectively. 
Testing  the  equation  for  pth  patch  on  mth  layer  using  /•„  „  \  («=1,  2,  ...,  XJ)  and 

yB^piy -ymp,v)  (v=l,  2,  ...,  V)  in  order,  will  generate  U(mp)+V(mp )  equations.  In  such  a  way, 

N  Km) 

testing  every  equation  for  the  corresponding  patch  results  in  ^  +  V(mp)\  equations  with 

m=4  p=  1 

unknowns  coefficients  of  the  same  number.  In  matrix  form,  we  have  AI  =  C  where  A  is  a  square 
matrix  (usually  called  impedance  matrix),  I  and  C  are  column  matrices  respectively  containing  the 
coefficients  of  the  current  expansion  functions  and  excitations.  A  typical  matrix  solver  like  LU 
algorithm  can  be  used  to  solve  for  I. 

Once  the  coefficients  of  the  current  expansion  functions  are  obtained,  all  patch  currents 
including  the  one  on  the  feed  can  be  calculated  using  (2).  A  de-embedding  scheme  similar  to  the 
one  in  [7]  can  then  be  used  to  obtain  the  voltage  reflection  coefficient  and  as  a  result  the  input 
impedance  Zjn.  The  total  radiated  far-fields  of  the  structure  can  also  be  obtained  by  integration  of 
the  current  distributions,  weighted  by  the  far-field  natural  modes,  over  the  surface  of  the  conducting 
patches.  For  the  sake  of  brevity  the  details  are  not  included  in  this  report. 

A.l.  Verification  of  the  Moment-Method  Code 

To  verify  our  developed  code,  we  first  simulate  a  single-layer  patch  antenna  using  our  own 
code  and  Ensemble.  The  patch  is  a  rectangular  one  operating  at  7.5  GHz  and  its  dimensions  are 
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15.81 14x12.7712  mm2.  The  substrate  thickness  is  0.79375  mm.  The  relative  dielectric  constant  and 
loss  tangent  of  the  substrate  are  2.2  and  0.001.  The  feed  line  is  2.46467-mm  wide  and  feeds  the 
patch  (along  a  radiating  edge)  at  a  distance  of  3.991  mm  from  a  non-radiating  edge.  Figure  10(a) 
and  (b)  show  respectively  the  real  and  imaginary  parts  of  the  normalized  input  impedance  in 
comparison.  Figure  1 1  compares  the  return  loss  data  from  the  two  codes.  It  is  seen  that  our  data  are 
in  good  agreement  with  those  from  Ensemble  except  for  a  frequency  shift  of  0.3%. 


Fig.  10:  The  real  and  imaginary  parts  of  the  normalized  input  impedance  of  a  single  patch  antenna. 


Figure  1 1 :  Return  loss  of  the  patch  antenna 

The  second  example  is  simulation  of  stacked  patches  in  a  two-layer  media.  This  structure  was 
modeled  and  measured  by  Dubost  and  Beauquet  [9].  The  first  patch  of  dimensions  6.4x6.4  mm2  is 
fed  at  the  center  of  an  edge  by  a  0.6-mm  wide  microstrip  line  and  printed  on  the  first  lossless 
substrate  with  er,=2.55  and  t,=0.8  mm.  The  second  patch  of  dimensions  5.6x5.6  mm2  is  printed  on 
the  second  substrate  with  £r2=2.55  and  t>=1.6  mm  and  arranged  symmetrically  with  its  sides  parallel 
to  those  of  the  first  patch.  The  real  and  imaginary  parts  of  the  normalized  input  impedance  of  the 
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antenna  structure  from  the  two  codes  are  shown  in  Figure  12(a)  and  (b)  for  comparison.  In  Figure 
13,  the  return  loss  data  from  the  codes  are  compared  with  the  measured  ones  from  literature.  Again, 
a  good  agreement  between  the  data  from  our  code  and  those  from  Ensemble  is  observed.  Also,  both 
sets  of  simulated  return  loss  data  exhibit  accordance  with  the  measured  ones.  This  stacked-patch 
antenna  is  actually  a  broadband  example.  Our  calculated  and  Dubost’s  measured  bandwidths 
(VSWR=2)  are  13.1%  and  13.5%,  respectively. 


Fig.  12:  The  real  and  imaginary  parts  of  the  normalized  input  impedance  of  a  stacked-patch  antenna. 


Figure  13:  Simulated  return  loss  in  comparison  with  measured  data  by  Dubost  and  Beauquet 


B.  Implementation  of  the  Evolutionary  Programming 

Evolutionary  programming  (EP)  [10]  is  well  suited  for  global  optimization  of  continuous, 
discrete,  and  mixed  parameter  problems  that  one  may  often  encounter  in  electromagnetics  and 
antenna  design.  Like  other  evolutionary  computational  techniques,  such  as  GA  [11]  and  ES  [12], 
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EP  is  a  heuristic  population-based  search  procedure  that  incorporates  random  variation  and 
selection.  However,  whereas  GA  uses  a  combination  of  crossover  and  mutation  with  the  former 
being  the  main  mechanism  of  change,  EP  uses  mutation  as  the  only  variation  operator.  The  mutation 
only  nature  of  EP  may  provide  a  powerful  tool  in  design  of  the  problem  specific  operators  and 
facilitate  easy  integration  with  available  apriori  knowledge  about  a  given  problem.  We  have 
experimented  with  two  different  mutation  operators,  namely  Gaussian  and  Cauchy,  in  the 
optimization  module  of  our  electromagnetic  optimization  engine.  In  this  part,  we  first  present  the 
details  of  two  modified  EP  algorithms:  meta-EP  with  Gaussian  mutation,  and  Fast  Evolutionary 
Programming  (PEP)  with  Cauchy  mutation.  These  algorithms  combine  the  mutation-based 
evolution  of  EP  with  the  self-adaptation  of  the  strategy  parameters  in  ES.  The  algorithms  are  then 
applied  to  the  challenging  constrained  optimization  problem  of  the  printed  Yagi-like  array  of 
stacked  dipoles. 

B.l.  Algorithms  for  meta-EP  and  FEP 

The  EP  algorithm  with  self-adaptive  mutation  operator  for  global  optimization  of  an  n- 

dimensional  objective  function  <|>(x ),  x  =  [*(]),*  (2) . *(»)]  consists  of  five  basic  steps: 

initialization,  fitness  evaluation,  mutation,  tournament  and  selection. 

Initialization:  An  initial  population  of  |i  individuals  is  formed  through  a  uniform  random  or  a 
biased  distribution.  Each  individual  is  taken  as  a  pair  of  real-valued  vectors,  Vi  \i} 

where  x,  =[^0)^(2) . *,  (»)]  and  rf,  are  the  n-dimensional  solution  and  its  corresponding  strategy 

parameter  (variance)  vectors,  respectively,  rj/s  are  initialized  based  on  the  user-specified  search 
domains,  e  {xmin ,  xn**}",  which  may  be  imposed  at  this  stage. 

Fitness  Evaluation:  A  fitness  function  F(xt)  is  assigned  to  each  vector  *.  in  the  population.  F(x)  in 
general  is  obtained  from  the  objective  function  <j>(3c)  by  possibly  imposing  some  random  alteration  p 
and  the  scaling  them  back  to  negative  real  values  by  a  scaling  function  a, 

F(x)=a(<|>(3c),  P)+P(x)  (3) 

where  P  is  a  collection  of  penalty  criteria  assigned  to  possible  optimization  constraints. 

Mutation:  In  EP,  mutation  is  the  only  variation  operator  used  for  generation  of  offsprings 
population  from  the  parents  population,  and  therefore  EP,  in  contrasts  to  Gas,  is  asexual  by  nature, 
i.e.,  each  single  parent  produces  one  single  offspring.  Design  of  efficient  mutation  operators  is 
presently  an  ongoing  topic  of  research  in  evolutionary  computation.  Here  we  present  two 
algorithms,  meat-EP  [12,13]  and  Fast  Evolutionary  Programming  (FEP)[14],  which  use  different 
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mutation  operators  in  the  evolution  process.  In  meta-EP,  each  parent  (*,,  q.)  creates  a  single 


offspring  (3c/,  q,')  by: 

Vu>  =  1+ fo/j  1 

(4) 

n,'o)  =  >1,0  1 

(5) 

for  j  =  0,1,2 . n,  where  x(j)  and  r\(j)  and  are  the  jth  components  of  the  solution  vector  and  the 

variance  vector,  respectively.  N(0,1)  denotes  a  one-dimensional  random  variable  with  a  Gaussian 
distribution  of  mean  zero  and  standard  deviation  one.  Nj(0,l)  indicates  that  the  random  variable  is 
generated  anew  for  each  value  of  j.  The  scale  factors  x  and  x’  are  commonly  set  to  (V*7J)  ‘  and 

(■&)'',  respectively,  where  n  is  the  dimension  of  the  search  space.  Self-adaptive  mechanism  of 

equation  (5),  borrowed  from  ES,  enables  the  meta-EP  to  evolve  its  own  variance  parameters  during 
the  search,  exploiting  an  implicit  link  between  internal  model  and  good  fitness  values.  The 
logarithmic  normally  distributed  process  for  the  variances  in  (5)  guarantees  positive  values  of 

standard  deviations.  The  global  factor  x’  N(0,1)  allows  for  an  overall  change  of  mutability  and 
guarantees  the  preservation  of  all  degrees  of  freedom,  whereas  the  factor  x  Nj(0,l)  allows  for 
individual  changes  of  the  variances  T|(j). 

In  FEP  algorithm  the  offsprings  are  still  generated  according  to  the  equations  (4)-(5),  but  with 
a  Cauchy  mutation  operator  replacing  the  Gaussian  mutation  in  (4),  i.e., 

'O' )  =  xt(J )  +  )  Cj ( 0 , 1 )  (6) 

where  C(0,1)  is  a  random  variable  with  a  Cauchy  distribution  operator,  G,  centered  at  the  origin  and 
with  the  scale  parameter  t  =  1 , 

Gt  0c)  =  —  +— tan-1(— )  ;  -°°<x  <°°  (7) 

2  it  t 

We  note  that  the  inverse  of  the  operator  in  (7),  needed  in  generation  of  the  random  variables  in  (6), 
is  given  in  closed  form.  It  has  been  shown  that  FEP  performs  much  better  than  EP  using  the 
Gaussian  mutations  for  many  unconstrained  optimizations  involving  multimodal  functions  with 
many  local  minima  as  well  as  in  nonlinear  signal  processing. 

Tournament:-  In  tournament  process  a  pairwise  comparison  with  respect  to  the  fitness  values  in  (3) 
over  the  union  of  parents  and  offspring  populations  is  conducted.  For  each  individual  sk  in  the 

union,  ke{l,....  2  p},  q  opponents  are  chosen  at  random  with  equal  probability  from  the  total 


18 


membership  2\i  of  the  union.  For  each  comparison,  if  the  individual  fitness  F  (ak)  is  no  greater  than 

the  opponent’s,  it  receives  a  “win”.  The  best  individual  is  guaranteed  a  maximum  ‘win’  score  of  q 
and  its  survival  to  the  next  generation.  We  note  that  this  tournament  process  differs  from  the  one  in 
Gas  and  is  ‘elitist’  in  nature. 

Selection:  The  JO.  individuals  out  of  the  union  of  (*,,  rj,ju(3c/,  Vi  e  {/,....  2  p}  with  the  most 

“win”  score  are  selected  to  be  the  parents  of  the  next  generation.  The  above  steps  in  the  fitness 
evaluation  through  the  selection  process  is  then  repeated  until  an  acceptable  solution  is  obtained  or 
the  prescribed  total  number  of  generations  is  reached. 

Figure  14  shows  the  flow  chart  of  our  electromagnetic  optimization  engine. 


Output 

Figure  14:  An  Electromagnetic  Optimization  Engine  for  Multi-layered  Printed  Antenna  Structures 
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C.  Examples  of  Gain  Optimization  of  the  Yagi-Like  Antennas 

To  optimize  the  gain,  the  Yagi-like  structure  is  represented  as  a  vector  consisting  of  the  length 
of  the  driver  element,  the  director  lengths,  the  relative  permittivity  and  the  dielectric  thickness  as 
given  below: 

X=[Ldri,  LSr{\) . LMr{N),  £r(l) . £r(N),  dl(l),...,dl(N)]T  '  (8) 

where  and  L^i),  er(i)  and  dl(i)  are  length  of  the  driver  element,  length  of  the  ith  director 
element,  dielectric  constant  and  thickness  of  the  ith  dielectric  layer,  respectively.  Length  of  the  vector 
in  (8),  in  view  of  the  mutations  (4)-(5),  is  in  general  n  =  6N+1.  For  the  gain  optimization  we 
construct  the  Fitness  function  in  (3)  with  a  =  -1  and  P  =  0, 

F(x  )  =  -  Gain®,  fc  3c)  +  Pm(x )  (9) 

m 

where  Gain  is  the  power  gain  in  (0,  <|>)  direction  obtained  from  the  moment  method  solution  of  a 

mixed-potential  integral  equation  for  the  current  distributions  on  the  strip  elements;  Pm  ,  m  =  1,2,.., 
are  the  penalty  criteria  for  violating  a  set  of  constraints  chosen  to  ensure  that  the  optimized  solutions 
are  practically  feasible. 

In  the  first  example,  we  present  the  performances  of  the  EP  and  FEP  algorithms  in  an 
unconstrained  optimization  of  a  five  layer  Yagi  array  with  air  dielectric  layers.  The  population  size 

and  the  number  of  opponents  were  set  to  |i  =  100  and  q  =  10,  respectively.  Optimization  was 

performed  with  respect  to  lengths  of  the  elements  and  thickness  of  the  layers.  Figure  15  shows  the 
mean  fitness-value  trajectory  of  the  best  population  member  at  each  generation  after  5  trial  runs 
were  performed.  As  seen  FEP  performs  better  than  meta-EP  in  terms  of  convergence  rate.  A  gain  of 
better  than  17dBi  was  obtained  after  155  generations. 
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Figure  15:  Optimization  of  a  5-layer  printed  Yagi  antenna  with  air  dielectric  using  two  mutation  operators 
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In  the  second  example,  a  Yagi-like  array  in  a  3  layer  isotropic  medium  was  optimized  with  (I  = 

20,  q  =  8  and  subject  to  the  constraints  of  D  =  dj+d2+d3  <  0.15X0  and  radiation  efficiency,  es  > 

95%.  The  optimization  parameters  in  this  case  were  lengths  of  the  elements,  thickness  of  the  layers 
and  the  dielectric  constants,  which  were  allowed  to  vary  in  the  range  of  2  to  4.  It  was  observed  that 

Eri  and  £r2  always  tended  to  2  while  tended  to  4;  this  is  consistent  with  the  high  gain  condition 

for  an  embedded  dipole  [2].  The  corresponding  mean  fitness  trajectories  for  meta-EP  and  FEP  after 
5  trials  are  plotted  in  Figure  16.  Meta-EP  performs  better  than  FEP  for  this  constrained 
optimization  case.  Gain  of  about  13  dBi  with  D  =  0.12  X0  and  es  =  99%  was  obtained  within  200 
generations.  It  is  noteworthy  that  the  same  structure  without  the  director  elements  has  a  gain  of  only 
about  7.5  dBi.  Gain  and  surface-wave  powers  of  this  structure  are  plotted  as  a  function  of 
normalized  frequency  in  Figure  17.  As  can  be  seen  at  the  maximum  gain,  the  power  coupled  into 
the  surface  waves  is  minimized. 
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Figure  17:  Gain  and  surface- wave  power  versus  normalized  frequency 
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IV.  DESIGN  OF  HIGH  GAIN  CIRCULARLY  POLARIZED  MICROSTRIP  PATCH 

ANTENNAS 

In  order  to  demonstrate  the  design  feasibility  of  multi-layer  high  gain  microstrip  antennas  for 
circular  polarization  (CP)  applications,  a  three-layer  rectangular  patch  antenna  is  designed  for 
MSAT  operation.  The  specifications  require  a  gain  of  larger  than  10  dBi  with  an  axial  ratio  (AR)  of 
less  than  1  dB  at  a  center  frequency  of  f= 1.643  GHz. 

A.  Design  Procedure 

We  have  used  a  nearly  square  patch  fed  by  a  coaxial  line  along  its  diagonal  as  shown  in 
Figure  18.  For  aDuroid  substrate  with  £r,  =  2.2  and  d,  =  1.59  mm,  the  dimensions  L  and  W  of  the 
patch  were  initially  selected  for  CP  according  to  [15]:  L  =  W  (1+1/Q,)  where  Q,  is  the  total  quality 
factor  of  patch.  Using  a  MOM-based  simulation  code,  the  feed  position  was  then  optimized,  in  the 
absence  of  any  superstate  layer,  until  a  low  VSWR  was  achieved.  Figure  18  shows  the  axial  ratio, 
gain  and  return-  loss  (SI  1)  of  the  initial  design.  At  the  frequency  of  1.643  GHz,  AR  =  0.278  dB 
and  VSWR  =  1.2.  The  realized  gain,  however,  is  only  about  7.1  dBi.  In  order  to  enhance  the  gain,  a 
high  permittivity  superstrate  layer  of  thickness  CI3  is  added  at  a  distance  of  d2  from  the  patch  surface 
as  shown  in  Figure  19.  A  high  gain  pattern  can  now  be  obtained  by  adjusting  the  thickness  d2  when 
d3  is  fixed  according  to  the  resonance  condition  [2]:  d}  =  \/4je^.  It  is  noteworthy  that  the 

presence  of  the  air  layer  allows  an  easy  mechanical  tuning  of  d2  in  a  practical  design  [16]. 
Resonance  frequency  of  the  original  patch  shifts,  however,  when  the  superstrate  layers  are  added 
which  hence  requires  adjustments  of  the  patch  dimensions  and  the  feed  position.  After  many 
simulation  trials,  we  have  finally  achieved  a  realized  gain  of  about  1 1.68  dBi  with  an  axial  ratio  of 
AR=1.08  dB  and  VSWR  of  1.25  at  f=1.643  GHz.  The  antenna  performance  as  a  function  of 
frequency  together  with  the  final  design  parameters  is  shown  in  Figure  19. 

B.  Performance  of  the  Designed  Antenna  in  a  Linear  Array  Environment 

It  may  be  desirable  to  also  form  an  array  of  such  high  gain  antenna  elements  for  even  higher 
gain,  low  side  lobe  and/or  steered-beam  applications.  Since  these  elements  have  relatively  thick 
dielectric  layers,  propagation  of  surface-waves  and  strong  mutual  couplings  may  adversely  affect 
the  gain,  axial  ratio  and  input  impedance  of  these  antennas.  To  investigate  the  coupling  effects  we 
have  analyzed  a  two-element  array  of  the  above  three-layer  design.  Figures  20-a  and  20-b  show  the 
gain  and  axial  ratio  of  the  array  for  the  E-plane  coupling  (i.e.,  when  the  elements’  W  sides  are 
facing  each  other)  and  the  H-plane  coupling  (i.e.,  when  the  elements’  L  sides  are  facing  each  other), 
respectively.  As  can  be  seen,  for  a  typical  0.5  \  spacing,  the  axial  ratio  has  degraded  from  1.08  dB 
for  an  element  in  isolation  to  more  than  5  dB  for  the  array.  Therefore,  for  these  multi-layered 
configurations,  one  needs  to  optimize  the  whole  array  for  both  high  gain  and  low  axial  ratio. 
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Figure  18:  Axial  ratio,  gain  and  return  loss  of  a  single  layer  CP  microstrip  patch  antenna 


-i 

vdi 


m**FmS?W**  .  .  t  -fa:-.*  -r:v: 

Axial  Ratio  ,  j 

™v„ _ _ | _ _  Gain  _ | _ —j#*?. 

**•  ■ — — 

|x  !  Jx  l 

*  1  Xi  yt  '  i 

i  j  ?  ; 


'  . . - ^  -J—  —  ■ 


|  j  ^  ? 


L=59.7  mm,  W=60.52  mm  !  X  1  : 

-»* - 4- — f^r- - — j-  Q  - -4* 

j  l  x  i  J? 

xn=39.6  mm,  dl =1.587 5  mm  — X~ — \ — X*+ — - — t -*t-- >- 

0  !  I  X..  i  X  < 

d.2-77 mm,  d3=24.5  mm  ^ 

Figure  19:  Axial  ratio,  gain  and  return  loss  of  a  3-layer  high  gain  CP  microstrip  patch  antenna 
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Figure  20:  Effects  of  mutual  coupling  on  the  gain  and  the  axial  ratio  of  the  array, 
(a):  E-plane  coupling;  (b):H-plane  coupling 


V.  GAIN  ENHANCEMENT  METHOD  FOR  MICROSTRIP  ANTENNAS  IN  UNIAXIAL 

ANISOTROPIC  MEDIA 


As  a  first  step  toward  the  extension  of  the  Yagi-like  concept  to  antennas  printed  in  multi-layer 
anisotropic  media,  we  have  analyzed  the  conditions  for  obtaining  high  gain  together  with  high 
radiation  efficiency  for  a  microstrip  dipole  antenna  embedded  in  a  multi-layer  anisotropic  medium. 

A.  Analysis  of  High  Gain  Conditions 

The  conditions  for  obtaining  high  gain  from  an  embedded  dipole  in  a  multi-layered  isotropic 
medium  have  been  discussed  in  [2, 3]  where  it  was  shown  that  it  is  possible  to  obtain  very  high  gain 
by  alternating  between  high  wave-impedance  (large  jir)  and  low-wave  impedance  (large  er )  layers. 
In  this  work  we  have  extended  the  above  gain  enhancement  method  to  the  case  of  a  microstrip 
dipole  printed  in  a  multi-layered  uniaxial  anisotropic  medium  with  permittivity  and  permeability 
dyadics  of  (Figure  21) 
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By  using  an  equivalent  cascaded  transmission  line  method,  explicit  expressions  for  the 
corresponding  far-fields  were  derived  and  the  conditions  of  obtaining  high  gain  together  with  high 
radiation  efficiency  (i.e.,  ratio  of  the  radiated  power  to  the  total  of  the  radiated  plus  the  surface-wave 
powers)  were  analyzed.  In  particular,  we  have  shown  that  by  a  proper  selection  of  anisotropy  ratios 
and  thickness  of  the  layers,  it  is  possible  to  obtain  two  beams  at  0E  and  0H  in  E  and  H  planes 
respectively.  More  explicitly,  for  the  case  of  a  three  layer  structure  in  Figure  21,  with 
F,  =  f2  ,  ^  =  ,  the  high  gain  conditions  are  derived  as: 


Figure  21:  A  Microstrip  Dipole  in  a  Multi-layered  Uniaxial  Medium 
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—  sin(0e)  = 

/Mn 

d^=d2=  1 

4«/l  . 

In  general,  the  condition  in  (1 1)  has  to  be  satisfied  for  each  uniaxial  layer  in  order  to  generate  two 
beams  exactly  in  0E  and  0H  directions.  One,  however,  can  still  obtain  two  beams  approximately  at 
0E  and  0H  if  only  the  first  two  layers  satisfy  (11).  In  that  case  the  third  layer  can  be  either 
anisotropic  or  isotropic.  More  specifically,  for  the  case  of  nonmagnetic  layers  the  high  gain 
conditions  can  be  reduced  to, 

4,2  =  -  ;  4=—^=  ;  V^sin(0£)  =  ^sin(0w)  (13) 


•sin  (6h) 


(ID 


sin2(0£) 

2  A 

;  4  =  7 

1  /■— \ 

CD 

ca 

a 

'Tsi 

I 

i _ 

4rt,3 

L  J 

(12) 


which  together  with  £0  »  1  produce  two  high  gain  beams.  Aside  from  a  normalization  factor,  the 
corresponding  far-fields,  approximated  near  the  main  beams  in  the  0E  and  0H  directions,  are  given 
by: 
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These  expressions  clearly  show  that  for  ea3  /  effl  »1,  Ee  and  E+  provide  two  beams  at  0E  and  0H  , 
respectively. 

The  above  formulation  is  valid  for  an  infinitesimal  printed  Hertzian  dipole.  In  order  to  analyze 
resonant  or  other  finite-length  printed  dipoles  in  multi-layered  uniaxial  media,  the  moment-method 
based  MPDE  code  was  extended  to  include  uniaxial  layers  by  modifying  the  corresponding  Green’s 
functions.  We  also  plan  to  extend  the  code  to  patch  antennas  of  arbitrary  shapes. 
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B.  Numerical  Results 

For  0E  =  60°  and  0H  =  30°  one  requires,  say,  8Z]=  8z2  =3  and  8^  =  8^=1.  Setting  the  third 
layer  to  be  isotropic  with  el3  =  ez3  =  26  ,  then  results  in  gains  of  about  17.1  dBi  in  H-plane  and 
about  9  dBi  in  E-plane  together  with  a  radiation  efficiency  of  80%  for  a  nearly  resonant  strip  dipole 
of  length  104  mm  at  f=  1  GHz..  The  input  impedance  is  =  154  -  j40  ohms.  The  corresponding 
patterns  are  shown  in  Figure  22-a.  The  minor  beam  pointing  error  in  the  E-plane  pattern  is  because 
the  equation  (11)  is  not  satisfied  for  the  third  layer.  The  gains  can  be  increased  by  increasing  the 
number  of  isotropic  layers  and  alternating  between  high  and  low  dielectric  constants  according  to 
the  resonant  conditions  in  [3,  17].  Figure  22-b  shows  the  patterns  for  a  five  layer  case  with  two 
anisotropic  layers  with  8zl=  8z2  =3  and  8tl  =  8^  =1,  and  three  isotropic  layers  with  e#  =26,  =2 

and  £,3  =26.  Gains  of  22.4  dBi  in  H-plane  and  17.1  dBi  in  E-plane  together  with  a  radiation 
efficiency  of  98%  are  obtained.  The  input  impedance  in  this  case  is  Zm  =  141  -  j46  ohms. 

Finally,  we  have  shown  that,  similar  to  the  isotropic  case  previously  reported  in  the  literature 
[18],  the  primary  radiation  mechanism  in  the  uniaxial  case  is  also  the  leaky-wave  excitation.  Details 
of  this  analysis  will  be  reported  in  the  next  annual  report. 


Angle  Angle 

(a)  (b) 


Figure  21:  Normalized  far-field  patterns  of  a  near  resonant  microstrip  dipole  antenna,  (a):  antenna  embedded  in  a  3- 
layer  uniaxial  medium;  (b):  antenna  embedded  in  a  5-layer  uniaxial  medium. 
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VI.  FUTURE  WORKS 


Our  research  works  on  the  multi-layered  printed  antennas  that  are  either  presently  in  progress 
or  are  planned  for  the  1999-2000  academic  year  include  the  followings: 

-  We  are  presently  in  the  midst  of  applying  our  electromagnetic  optimization  engine  for  gain 
optimization  of  microstrip  line-fed  Yagi-like  patch  antennas  and  arrays  for  linear  as  well  as  circular 
polarization.  Other  configurations  planned  for  optimizations  include  proximity-coupled-fed  and 
coaxial-fed  microstrip  patch  antennas.  We  also  plan  to  fabricate  and  experimentally  verify  selected 
optimized  designs  for  such  applications  as  in  GPS  and/or  MS  AT. 

-  We  plan  to  use  our  electromagnetic  optimization  engine  to  develop  novel  schemes  for  shape 
optimization  of  printed  antennas  for  various  multi-function  applications  such  as  multi-band,  wide¬ 
band  and/or  multi-polarization  operations.  We  also  plan  to  fabricate  and  measure  a  few  shaped- 
optimized  microstrip  antennas  in  order  to  validate  any  novel  design  concepts,  which  may  result  from 
this  investigation. 

-We  plan  to  extend  our  electromagnetic  optimization  engine  to  include  Yagi-like  stacked  patch 
antennas  printed  in  a  multi-layer  anisotropic  medium.  The  optimization  engine  will  then  be  used  to 
investigate  these  antennas  for  high  gain  multi-beam  applications. 

-  We  plan  to  analyze  the  Yagi-like  antennas  in  terms  of  the  leaky-wave  excitation  in  order  to  better 
understand  the  physics  of  their  high  gain  radiation.  The  analysis  will  be  performed  for  both 
isotropic  as  well  as  uniaxial  anisotropic  cases. 
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